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COEXISTENCE IN STOCHASTIC SPATIAL MODELS^ 

By Rick Durrett 

Cornell University 

In this paper I will review twenty years of work on the question: 
When is there coexistence in stochastic spatial models? The answer, 
announced in Durrett and Levin [Theor, Pop. Biol. 46 (1994) 363- 
394], and that we explain in this paper is that this can be determined 
by examining the mean-field ODE. There are a number of rigorous 
results in support of this picture, but we will state nine challenging 
and important open problems, most of which date from the 1990's. 

Introduction. There is an incredible diversity of species that coexist in 
the world. At the Botanic Garden in Singapore one can see 1000 species of 
orchids. These are cultivated, of course, but if one examines the food web in 
a small lake one finds dozens of species coexisting. An important problem in 
ecology is to identify mechanisms that permit the coexistence of species. In 
this paper we will examine that question in the context of stochastic spatial 
models. In these interacting particle systems, space is represented by the 
d-dimensional integer lattice Z*^. With ecological problems in mind, we will 
usually take d = 2. 

Our story begins with a very simple and natural model, but one that is 
still not well understood. 

Example 0. Competing contact processes. 

• Each site in 1? can be in state = vacant, or in state i = 1, 2 to indicate 
that it is occupied by one individual of type i. 
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Fig. 1. Competing contact process ODE. 

• Individuals of type i die at rate 6i and give birth at rate Pi. Here, at rate 
A means that these events happen at times of a rate A Poisson process. 

• A type i born at x goes to x + y with probabihty Pi{y). If x + y is vacant 
it changes to state i, otherwise nothing happens. 

When there is only one type this reduces to the system introduced by Har- 
ris (1974). After several decades of work, this model is very well understood. 
See Liggett (1999) for a survey. 

If we assume that the states of adjacent sites are independent then the 
fraction of sites Ui in state i = l,2 satisfies 

^ = Piui{l -ui- U2) - 5iui, 

(1) 2, 

—— = I32U2{1 -Ui- U2) - S2U2. 

at 

This is called the mean-field ODE, because if we consider the system on N 
sites with a uniform dispersal distribution then in the limit as — > 00 the 
densities converge to this limit. In the spatial model adjacent sites are not 
independent. However, writing and analyzing the mean-field ODE is a good 
first step in guessing what the system will do. 

In (1) dui/dt = when {1 — ui — U2) = Si/Pi. These lines are parallel, so 
they either do not intersect or coincide. Figure 1 shows the mean-field ODE 
when Pi = 4, /32 = 2, and 5i = 62 = 1- In this case all solutions starting from 
a point with ni > converge to (3/4,0). 

Neuhauser (1992) proved the following result: 

Theorem 1. // the dispersal distributions pi =P2=P, death rates 5i = 
82 = 5, and birth rates Pi > P2 then species 1 out competes species 2. That is, 
if the initial condition is translation invariant and has P{^o{x) = 1) > then 
P(6(x)=2)^0. 
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Fig. 2. Simulation of competing contact process. /3i = 3.9, Si = 2 (black) versus [32 — 2, 
Si = 1 (gray). The picture is a snapshot of part of the grid at time 300. 

Ideas behind the proof. We construct the process from a "graphical 
representation." For each site x there is a rate 5 Poisson process D^, n > 1 
that kihs the particle at x (if there is one). For each x and y, there are Poisson 
processes B^^y , n>l and A^'^, n > 1 with rates j32p{y) and — (32)piy)- 
The first causes births from x to x + y if x is occupied and x + y is vacant. 
The second causes births x to x + y if x is occupied by a 1 and x + y is 
vacant. Using this construction and working backwards in time, Neuhauser 
(1992) was able to show that the extra arrows reserved only for the I's meant 
that if the site was occupied then it would be a 1 with high probability. We 
refer the reader to Neuhauser (1992) for the details, which are somewhat 
complicated. Our reason for giving this sketch is to make clear that the 
proof requires 5i = ^2 • D 

In the mean-field ODE, only the ratios l3i/6i matters, so it is natural to 
guess. 

Problem 1. Show that species 1 outcompetes species 2 holds if the 
dispersal distributions are the same and f3i/5i > (32/ The simulation in 
Figure 2 gives some support for this conjecture. 

The behavior in Problem 1 is what biologists expect based on the Com- 
petitive Exclusion Principle. A version of this can be found in work of Levin 
(1970). Consider an ODE of the form: 



= Uifi{zi,...,Zm), l<i<n. 
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Fig. 3. Host pathogen mean-field ODE. 

Here the Zi are resources. In the previous model there is one resource: zi = 
1 — ui — U2 free space. A special case is shown in Figure 3. Omitting the 
details of the assumptions on the fi in the mathematical result, the principle 
as biologists use it is: 

Theorem 2. If n > m then no stable equilibrium in which all n species 
are present is possible. 

In words, the number of coexisting species is smaller than the number of 
resources. 

To try to build some excitement about Problem 1, we note that Chan 
and Durrett (2006) proved that a "fugitive species" that disperses at a fast 
rate, according to a truncated power law distribution, can coexist with a 
superior competitor with a nearest neighbor dispersal distribution. This does 
not contradict the competitive exclusion principle, because in addition to 
single site deaths their model has forest fires, which destroy large squares of 
occupied sites at a small rate. Thus the model has a second type of space, 
"recently disturbed space," and is entitled to have two coexisting species. 

As announced in the abstract, one goal of this paper is to explain the 
idea of Durrett and Levin (1994) that one can determine whether coexis- 
tence happens in the stochastic spatial model by examining properties of 
the mean-field ODE. The discussion is divided into three cases according to 
the properties of the ODE. 

Case 1. Attracting fixed point. When the mean-field ODE has an at- 
tracting fixed point with all components positive, we expect coexistence in 
the spatial model, that is, there is a stationary distribution which concen- 
trates on configurations that have infinitely many sites occupied by each 
species. 
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Example 1.1. Grass hushes trees. In this variant of the contact process 
there is a hierarchy of types. In hindsight this is a very natural model. 
However, it owes its invention to talking to Simon Levin about successional 
sequences in a forest. 

• Each site in 1? can be in state = grass, 1 = bush, 2 = tree. 

• Individuals of type i die at rate 6i and give birth at rate 

• A type i born at x goes to x + y with probability Pi{y). If a; + y is in state 
j < i it changes to state i, otherwise nothing happens. 

The mean field ODE is 

-ui- U2) - 5iui - P2U2U1, 

(32U2{1 - Ui) - 52U2. 

If P2> S2, the equilibrium frequency of 2's, U2 = {P2 — ^2)1 1^2- Given this, 
one can solve for u\ and see when it is positive. However, it is better to 
approach the question by examining when the "I's can invade 2's in equi- 
librium." The phrase in quotes means that if the 2's are in equilibrium and 
the I's are at a small density then the density of I's will increase. Ignoring 
the possibility that a 1 will see another 1 nearby, the condition is: 

(3) /?! • -5- > (^1 + /32 • — 5 • 

P2 P2 

The left-hand side gives the rate at which I's give birth onto O's, while on 
the right, the first term is the death rate of I's and the second is the rate at 
which they are eliminated by births of 2's. It is easy to check that if > ^2 
and (3) is satisfied then (2) has an equilibrium will all components positive. 

When 5i = = 1, (3) becomes /3i > Z?! > 1- For simplicity, we will con- 
sider only this case. Durrett and Swindle (1991) have shown 

Theorem 3. > /Jf > 1 th,en when pi is uniform on{x:0 < \\x\\ < L} 

and L is large, there is a stationary distribution ^12 that concentrates on 
configurations with infinitely many I's and 2's. 

Our hammer: the block construction. The survival of the 2's is 
not a problem because they are a contact process and don't feel the presence 
of the I's. To prove that the I's can persist in the space that remains, we 
use a "block construction," which consists of comparing the particle system 
with a mildly dependent oriented percolation in which sites are open with 
probability close to 1. 

For an account of this method see my St. Flour lecture notes, Durrett 
(1995), or return to the first application, Bramson and Durrett (1988), for a 



(2) 



dui 
~dt' 

du2 
"dt 
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very simple example. At an intuitive level, what one shows is that "one pile 
will make two piles with high probability." 

Here, pile is an undefined term and the phrase is short hand for the main 
idea behind the construction. To be more precise, in d = 1 one might show 
that L and T can be chosen so that if [— L, L] is "good" at time 0, then 
[— 3L, —L] and [L, 3L] will be good at time T with probability 1 — e. In this 
case, good might be there are not too many 2's in [— 2L, 2L] and there are 
enough I's in [—L, L\. To get a finite range of dependence between the events 
in the construction, we need to estimate the probability of the good events 
if we assume that all sites in [—kL^ kVf are always occupied by 2's. 

If e < Efc, that is, a constant that depends on the range of dependence, 
then facts about oriented percolation give a lower bound on the density of 
sites occupied by I's. Taking the Cesaro average of the distribution from 
time to time i, and finding a convergent subsequence produces the desired 
stationary distribution. For more details on the last point, see Liggett (1985). 
□ 



The block construction technology has improved quite a bit since 1991, 
so at this point it should be fairly routine to do the following: 

Exercise. Show that if (5i = ^2 = 1, /32 > 1 and /3i < 02 then the I's die 
out when the range is large. 

In the setting of Theorem 3, in addition to existence of a stationary dis- 
tribution, we have a uniqueness result proved by Durrett and Moller (1991). 

Theorem 4. Suppose 5i = ^2 = 1, /?i > /?! > 1- ^/ ^he range is large then 
whenever the I's and the 2's do not die out then the process converges to 

^12- 

There are stationary distributions fii and fi2 with only I's and 2's respec- 
tively. By results for the one-type contact process, these are unique if one 
specifies that there is no mass on the all O's state. In addition there is the 
trivial stationary distribution /iq that assigns mass 1 to all O's. The con- 
vergence result in Theorem 4 when combined with results for the one-type 
contact process implies that all stationary distributions are convex combi- 
nations of Ho, fii, fi2 and /_fi2. 
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Interlude: fast stirring. There are, at this point, a number of coexistence 
results for particle systems with long range interactions: Durrett (1992), 
Durrett and Schinazi (1993), Durrett and Neuhauser (1997), Durrett and 
Lanchier (2008), etc. However, the proofs of these results are done on a case 
by case basis. Things are simpler if, instead of long range, we assume that 
there is "fast stirring": suppose that for each pair of nearest neighbors x 
and y, at rate exchange the values Ct{x) and ^t(y). In this case there is 
a general result. 

Theorem 5. Suppose there is a function (j) that (i) decreases along 
solutions of the mean-field ODE, (ii) is convex, and (iii) (p{u) — > oo when 
mini Ui —>■ 0. Then there is coexistence in the model with fast stirring when 
e<eo. 

Durrett (2002) applies this result to a wide variety of systems: epidemics, 
predator-prey models, predator mediated coexistence, etc. 

Sketch of the proof of Theorem 5. Suppose that the mean-field 
ODE is 

du 

Conditions (i)~(iii) and a few lines of calculus implies, see pages 102-103 
of Durrett (2002) that if h = (p(u) and (p is the time derivative of cp along 
solutions of the ODE then 

dh . , , 

— < Ah + (p(u). 
dt 

Since (j){u) < 0, this is Brownian motion with killing. Using this we can show 
that there are constants 6,c> and T < oo, so given starting conditions 
with Ui{0,x) > 77 > for X G [—L,L], solutions of the "mean-field PDE" 

— = An + /(n) 

have miuj Ui{t, x) >5 for t>T, \x\ < ct. As e —>■ the particle system on eZ'^ 
converges to the solution of the PDE. Using this with the result for the PDE, 
we have shown that "one pile will make two piles with high probability" and 
the result follows from the block construction. □ 



Problem 2. Formulate and prove a similar general result for systems 
with long-range interactions. 
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Example 1.2. Host-pathogen models. In Durrett (2002) it was shown that 
predation can cause two competing species to coexist. Durrett and Lanchier 
(2008) have shown that coexistence can occur if there is a pathogen in one 
species. In the next models 1 and 3 are the two species, while 2 is species 
1 in the presence of a pathogen. Letting /j be the fraction of neighbors in 
state i, the rates are 

1^2 a/2 
2^1 72(/l + /2) 
3^1 73(/l + /2) 

1 ^ 3 71/3 

2 ^ 3 72/3 

The first two rates say that if there are no 3's and we think of 1 = vacant 
and 2 = occupied then the I's and 2's are a contact process. To explain the 
last four rates: at rate 7^ individuals of type i are replaced by the offspring 
of a randomly chosen neighbor. If the neighbor is type 1 or type 3 then the 
offspring has the same type as the parent. However, if the neighbor is type 2, 
the offspring is of type 1 because the pathogen is not passed into the seeds. 

The mean-field ODE is 

^ = {Ui + U2)ij2U2 + 731^3) - aUiU2 - JlUiUs, 

(4) —T-=aUiU2-j2U2, 

at 

^ = 1*3(71^^1 + I2U2) - Isu-siui + U2). 
As we will explain in a moment this leads easily to the following 

Theorem 6. Suppose 71 < 73 < 72 < a and 

(5) 71— +72( 1 - — ) > 73 

a \ a / 

then there is coexistence for large range. 

Keys to the proof. Here the I's and 2's are a contact process so on 
the boundary U3 = 0, = 72/a is an attracting fixed point. The displayed 
condition says that the 3's can invade the I's and 2's in equilibrium. The 
key to the proof is using an understanding of the ODE to show that if the 
density of some type becomes small then a sequence of events will occur that 
results in all of the densities being bigger than some e. The two observations 
which lead to this are: (i) the boundary ui = is not a problem since 2's 
give birth to I's; (ii) on the boundary U2 = 0, I's and 3's are a biased voter 
model in which I's outcompete 3's. □ 
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Figure 4 gives two simulations. It is clear that one cannot have coexistence 
if 7i and 72 are both > 73 or both <73. The next problem address the 
remaining case: 

Problem 3. Show that coexistence is not possible in the host-pathogen 
model if 72 < 73 < 7i • 

In this case the pathogen is called a mutualist, since it decreases the rate 
at which the species is replaced. To explain why the result in Problem 3 
should be true, note that if we start with the I's and 2's in equilibrium and 
a small density of 3's then once the invasion of the 3's starts the fraction of 
2's gets smaller, and the 3's have an even bigger advantage over the I's and 
2's. It is not hard to check in this case that there is no interior fixed point. 
For more results and problems about host-pathogen models see Lanchier 
and Neuhauser (2006). 

Case 2. Two locally attracting fixed points. In this case, the limiting 
behavior of the ODE depends on the initial density. However, this is not the 
expected behavior for the particle system, and in this case the outcome of 
competition is dictated by the behavior of the PDE. The reason is that even 
if the initial distribution is translation invariant and hence has a well-defined 
density, there will be regions of space where the density of I's is close to 1 
and others where it is close to 0. To explain what we expect to happen in 
the particle system, we consider an example. 




Fig. 4. Host pathogen systems with coexistence and no coexistence. 1 = black, 2= white, 
3 = gray. 
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Example 2.1. Sexual reproduction. The flip rates are as follows: 

• 1 ^ at rate 1, 

• ^ 1 at rate f3k{k — \)/n{n — 1) if A; of the n neighboring sites are occu- 
pied. 

In words, at rate /3, a vacant site picks two of its neighbors at random 
and become occupied if they both are. 
The mean-field equation is: 

(6) '^ = -u + l3u^{l-u)=u{-l + l3u{l-u)). 

Remembering that u{\ — u) is maximized at 1/2, where the value is 1/4, we 
see that there are nontrivial fixed points pi < p2 if and only if /3 > 4, while 
if /3 = 4, 1/2 is a double root. 

At this point one might think that in the presence of fast stirring, the 
critical value for survival of the process /3c ~ 4 but that is not correct. To 
determine the asymptotics of the critical value you have to consider the 
mean-field PDE: 

(7) — = An + <7(n), 

where g{u) = u{—l + /3ti(l — u)). 

A solution of (7) of the form u{t,x) = w{x — ct) with ^/{—oo) = p2 and 
w{+oo) = is called a traveling wave. In order to be a solution w must 
satisfy 

—cw' = w" + g{w). 
Multiplying by w' and integrating from — oo to oo 

{w')'^dx= / w"w'dx+ / g{w)w' dx. 



The first antiderivative on the right is (?i/)^/2 which vanishes at — oo and 
oo. Changing variables y = w{x) in the second we have 

V) dx= g{y)dy. 
Jo 

Thus the sign of c = the sign of Jq^ g{y) dy. When /? = 4.5, the cubic has 
roots at 2/3, 1/3, and and symmetry dictates that J^^ g{y)dy = 0. Thus 
c > if and only if /3 > 4.5. 

Durrett and Neuhauser (1994) have shown: 

Theorem 7. // we introduce fast stirring at rate then in the sexual 
reproduction model, ^-5 as e — > 0. 
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Sketch of the proof. The key is the PDE fact that if > and the 
initial condition is u{0, x) > pi + rj for |x| < L and L> L^^ then for any 5 > 
we have u{t, x) ~ p2 for |x| < {c — 5)t where c is the wave speed defined above. 
Combining this with the convergence of the particle system with fast stirring 
to the PDE, we have the "one pile makes two piles with high probability" 
needed for the block construction. □ 

There are many situations in which we can prove the existence of station- 
ary distributions, but convergence results like the one for grass-bushes-trees 
in Theorem 4 are rare. One exception is the result for multicolor contact 
processes in Durrett (1992). One would hope to develop general methods 
for proving uniqueness, but it seems sensible to start with a concrete case. 

Problem 4. Consider the sexual reproduction model with (3 > 4.5. 
Show that when the stirring rate is large there is a unique nontrivial 
stationary distribution. 

Example 2.2. Catalyst. Moving away from ecology, our next system is 
a model for the catalytic converter in a car's exhaust system. States are 
= vacant, 1 = CO (carbon monoxide attached to the surface), 2 = (single 
oxygen atom attached to the surface). The rates are as follows: 

• ^ 1 at rate p, 

• a pair of neighboring O's — > 22 at rate q/A, 

• adjacent 12 00 at rate r/4 (reaction to form C02). 

In this model all I's and all 2's are absorbing states corresponding to 
poisoning of the catalyst surface. In order for the catalytic converter to 
work and turn CO into CO2 there must be coexistence in the spatial model. 
Ziff, Gulari and Barshad (1986) considered the case in which r = 00 and 
q/2 = 1 — p (the latter condition can be imposed by scaling time). Their 
simulations shows coexistence for 0.389 <p< 0.525. 

Since I's land at rate p and two 2's land at rate <q = 2(1 — p), it seems 
clear that the system converges to all I's when p > q. There is a simple 
argument, see Theorem 1 in Durrett and Swindle (1994), which shows that 

p > q then P(^t(x) = 0) — > and if x and y are neighbors, P(^j(3;) = 
1,6(2/) = 2)^0, so 

P(6(x) = 1 on [-K, Kf) + P(6(x) = 2 on [-K, Kf) ^ 1 

but we do not know how to prove that if we start from the = configuration, 
the system converges to all I's. It is not hard to show, see Theorem 2 in 
Durrett and Swindle (1994), that the system converges to all 2's for small 
p. However, it is much more interesting to 
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Problem 5. Find < pi < p2 < 1 so that coexistence occurs in the 
catalyst model for p G (pi,P2)- 

Simulations suggest that the density of O atoms in equilibrium drops to 
discontinuously at the upper critical value, but proving this is a much harder 
problem. 

Bramson and Neuhauser (1992) have proved coexistence when the 02's are 
replaced by an x polymer, which needs a vacant N x N square to land, 
and is sufficiently large. Durrett and Swindle (1994) proved coexistence 
in the original model by introducing fast stirring. The mean-field PDE is: 

Aui + p{l - Ml - U2) - ruiU2, 

Au2 + q{l - ui - M2)^ - ruiU2. 

p < q, the ODE has four fixed points: two stable (1, 0) and (a, /?) and two 
unstable: (0,1) and (/?, a), where 

^0 {q-p)^ ^/[q-p)'^ - Aqp^/r 

a<p= . 

2q 

See Figure 5 for an example. 

The PDE results that were routine for the sexual reproduction model are 
now difficult. To prove the existence of a traveling wave with u(— 00) = 
(q,/3) and n(oo) = (1,0) one goes to the four dimensional phase plane: 
{ui,Ui,U2,U2), and looks for a curve connecting (a, 0,(3,0) and (1,0,0,0) 
which will exist only for one value of the speed c. Fortunately this was done 
previously by Volpert and Volpert (1988). With the existence of a traveling 
wave established the next step is to prove a convergence theorem for the 



(8) 



dui 

du2 
'dT 
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Fig. 6. Colicin mean field ODE. 



PDE, which can be done with comparison techniques because of monotonic- 
ity properties of system {ui,—U2). Once the PDE result is estabhshed the 
rest is a routine apphcation of the block construction. 

Example 2.3. Colicin. Durrett and Levin (1997) considered a competi- 
tion between two types of E. coli, one of which produces colicin (a chemical 
that kills other E. coli): 



Birth 


Rate 


Death 


Rate 


0^1 


Ml 


1^0 




0^2 


M2 


2^0 


^2 + 7/1 



Here the rates are like the two-species contact process, except for the 
7/1 in the death rate 2 — > 0, which comes from I's killing 2's with colicin. 
For simplicity we suppose that the basic death rates are equal 61 = 62 = ^■ 
Having done this it is natural to suppose that /3i < /32, for otherwise it is 
clear that the I's will out compete the 2's. 

The mean-field ODE is: 

^ = (3iui{l -Ul- U2) - 5iui, 

(9) 2,, 

—— = (32U2{1 - Ul- U2) - U2{52 + lUi). 

at 

A little algebra shows that when 

X / A <^2 ^ ^1 -^2 + 7 

di<(ii and -5- < -5- < ^ , 

P2 Pi P2 + 7 

the mean-field ODE has an interior fixed point but it is unstable. Figure 6 
shows the situation when f5i = 3, 71 = 2.5 and /?2 = 4. 

Figure 7 gives the density versus time in the system on a 200 x 200 grid. 
To emphasize that the behavior is different from the ODE, we start the gray 
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colicin producer (I's) at a small density. By time 1000 it has eliminated the 
black colicin sensitive strain (2's). The other panel shows the state at time 
600. Note that the two types have segregated. A movie would show that the 
interface moves in a direction that favors the I's. 

Problem 6. Show that coexistence is not possible in the colicin model 
when /3i < P2 and 61 = 62 = 1- 

Case 3. Cyclic systems, periodic orbits. In this case, we see coexistence 
with significant spatial structure. The pictures are pretty but the problems 
are hard. 

Example 3.1. Multitype biased voter model. Each site can be in state 1, 
2, . . . ,k, and j — > i at rate fiXij. In words, i's eat j's at rate Xij. The mean 
field ODE is 

j 

Silvertown et al. (1992) who were interested in the competition of grass 
species, studied the five species case in which 

0.09 0.32 0.23 



A,: 



0.08 
0.06 
0.02 
V0.O2 




0.06 
0.06 
0.03 



0.16 


0.05 
0.05 



0.06 
0.44 


0.03 



0.37\ 
0.09 
0.11 
0.03 
/ 



This example is not very interesting because Xij > Xji for 2 < j < 5, so using 
ideas of Grannan and Swindle (1990) and their improvement by Mountford 





Fig. 7. Simulation of colicin model. Producers (I's) are gray; sensitive strain (2's) is 
black. 
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and Sudbury (1992) one can show that if is the event I's are ahve at 
time t then P{Al ,^,t{x) 7^ 1) ^ so if the I's do not die out they take over 
the system. The key to the proof is that if 9 is small 



Durrett and Levin (1998) studied the cychc case in which A13 = 
A21 = /?2 5 A32 = Ps, and the other Xij = 0. This system with Pi = 1 and the 
corresponding discrete time deterministic cellular automata had been stud- 
ied earlier by Bramson and Griffeath (1989), Fisch, Gravner and Griffeath 
(1991) and Durrett and Griffeath (1993). The mean-field ODE has equilib- 
rium: Pi = /3j_i/(/3i + (32 + Ps) where z — 1 is computed modulo 3 with the 
result in {1,2,3}. Around this fixed point are concentric periodic orbits. See 
Figure 8 for an example. To prove mathematically that this occurs, write 
H{u) = XliPjlog^i check that H{u) is constant along solutions of the 



Problem 7. Show that coexistence always occurs in the cyclic case of 
Silvertown's model. 

For partial credit show that coexistence can occur for some parameters. 
Figure 9 gives a proof by simulation. 

Interlude: lizard love. Systems with a cyclic relationship exist in na- 
ture. In the side-blotched lizard [Uta stansburiana) , males have one of three 
throat colors, each one declaring a particular strategy. Dominant, orange- 
throated males establish large territories within which live several females. 
But these territories are vulnerable to infiltration by males with yellow- 
striped throats — known as sneakers — who mimic the markings and behav- 
ior of receptive females. The orange males can't successfully defend all their 



is a submartingale. 



z-4tiz)=l 



ODE. 




Fig. 8. Cyclic particle system with /3i — 0.3, /32 = 0.7, f3a = 1.0. 




Fig. 9. Cyclic particle system with /3i = 0.3, — 0.7, 13-^ — 1.0. 1 's are black, 2's white, 
3's gray. 



females against these disguised interlopers, who cluster on the fringes of the 
territories held by the orange lizards. 

However, a large population of sneakers, which have no territory of their 
own to defend, can be quickly overrun by blue-throated males, who defend 
territories large enough to hold just one female. Sneakers have no chance 
against a vigilant, blue-throated guard. But once the sneakers become rare, 
powerful orange males flourish, grabbing territory and females from the blue 
lizards. Now, the blue males lose out. See Sinervo and Lively (1996) for more 
on these lizards. 

Example 3.2. Three species colicin. Durrett and Levin (1997) considered 
an E. coli competition model with rates 



Birth 


Rate 


Death 


Rate 


0^1 


A/i 


1^0 




0^2 


hh 


2^0 




0^3 




3^0 


h +7i/i + 72/2 



Here, I's and 2's are colicin producers, while 3 is colicin sensitive. In the 
two species system (we conjecture) there is no coexistence, but as we will 
see coexistence is possible with three species. 

Consider for concreteness, the situation when 5j = 1, /3i = 3, /32 = 3.2, 
[3^ = 4.0, 71 = 3 and 72 = 0.5. In this case the 2's beat the I's since they have 
a larger birth rate, the 3's beat the 2's since the colicin they make is not 
nasty enough, while the I's beat the 3's. Thus again the three competitors 
have the same relationship as in the child's game paper-rock-scissors. 

The mean-field ODE is similar to (9): 

^ = /3itii(l -U1-U2- U3) - 6iui, 
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(10) = /32U2{1 - Ui - U2 - Us) - 62U2, 

^ = PsUsil -U1-U2- M3) - Uz{5s + 71M1 + 72ti2)- 

Figure 10 gives a picture of the mean-field ODE in the concrete case 
considered above as we look down into the tetrahedron Ui > 0, ui + U2 + us < 
1. On the (ui, 0, u^) and (0, ^2,^3) faces we see the colicin ODE, while on the 
(?xi,«2,0) face we have the competing contact process. As in the competing 
contact process comparing the first two equations shows that there is fixed 
point with U1U2 > when Pi/Si^ (32/^2- 

Figure 11 gives a simulation. The graph gives the numbers of the different 
types on a 200 x 200 grid, while the picture gives a snapshot of part of the 
system at the final time. 

Problem 8. Show that coexistence can occur in the three species colicin 
model. 



•1 




Fig. 1L Three species colicin simulation: I's black, 2's darker gray, 3's lighter gray. 
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Coexistence has been verified experimentally by Kirkup and Riley (2004). 
They began with a sensitive strain (S) of E. coli, introduced colicin plasmids 
to make a colicin-producing strain (C), and exposed the sensitive strain to 
the colicin-producer to obtain a resistant strain (R). A number of unlucky 
mice were then chosen to have the competition drama play out in their 
colons. After reading the paper, I think I will stick to computer simulations. 
In four cases, the experiment had to be discontinued because the mice were 
fighting or several mice died. 

Example 3.3. Spatial Prisoner's Dilemma. This time we allow multiple 
hawks T]t{x) and doves Ctix) at each site. 

• Migration. Each individual at rate migrates to a nearest neighbor. 

• Death due to crowding. Each individual at x dies at rate K{r]t{x) + Q{x)). 

• Game step. Let pt{x) be the fraction of hawks in the 5x5 square centered 
at X. Hawks give birth (or death) at rate apt{x) + 6(1 — pt{x)), doves at 
rate cpt{x) + d{l-pt{x)). 

Here, "birth (or death)" means that if the quantity is positive it is a birth 
rate, but if it is negative it is —1 times a death rate. 

An interesting choice for the game matrix is a Prisoner's Dilemma 





H 




D 


H 


a = -0.6 


b 


= 0.9 


D 


c= -0.9 


d 


= 0.7 



The H strategy dominates D, so it is the better choice, but the payoff 
for (D,D) is better than that for {H,H). This is the Prisoner's Dilemma 
"paradox." If everyone played D then the world would be a nice place, but 
this leads to the temptation to play H and increase your payoff. 




Fig. 12. Hawks-Doves ODE. 
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In a homogeneously mixing population the densities of Hawks (u) and 
Doves (v) would evolve according to 



In the concrete example under consideration, we have the following behavior. 
On the vertical axis one can see that in the absence of Hawks, Doves reach 
an equilibrium. However, when both are present, the ratio of Hawks to Doves 
increases until the population crashes to 0. See Figure 12 for an example. 

Simulations suggest that in the spatial model Hawks and Doves can coex- 
ist. Intuitively this occurs because the dynamics of the ODE happen locally, 
but when the Hawk population crashes to then there are some Doves left 
behind to recolonize space, but when the Doves grow to a significant density 
then the remaining Hawks again take over. 

Problem 9. Show that coexistence can occur in the spatial Prisoner's 
Dilemma. 

Credits. Since I first investigated the right edge of the one dimensional 
contact process in 1980, I have written 92 papers (out of my current total of 
170) on the topic of stochastic spatial model. On 21 papers, I am the only 
author, but most of those are survey articles and conference proceedings, so 
much of the credit for my work should go to my collaborators. Listed in order 
of their multiplicity they are: (12) Ted Cox; (10) Simon Levin and Linda 
Buttel; (6) Maury Bramson; (4) David Griffeath, Roberto Schonmann; (3) 
Tom Liggett, Ed Perkins, Rinaldo Schinazi, Jeff Steif; (2) Wan-ding Ding, 
Nicolas Lanchier, Xiu-fang Liu, Iljana Zahle; (1) Ben Chan, Larry Gray, 
Paul Jung, Anne Moller, and Mateo Restrepo. In addition I would like to 
thank two referees, and my current postdoc John Mayberry for helping me 
to improve the paper. 
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